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ABSTRACT 

Scattering  problems  are  considered,  from  the  point  of  view 
of  the  mathematical  questions  involved,   isach  as  the  boundary  condi- 
tions and  uniqueness   of  the  solutions,  and  the  determination  of  a 
Greenes  function.     For  sinrplicity,  the  treatment  Is  confined  to  one- 
particle  one-dimensional  problems. 
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I,   Introduction 

A  survey     of  the  voluminous  literature  on  scattering  suffices   to   show 
that  one-particle  scattering  problems,    or   those  which  are  effectively  one-particle, 
e.g.    two-particle  scattering  in  the  center-of-aase   system,    are  apparently  well  for- 
mulated and  the  theoiy    is  in  agreement  with  experiraent.      In  many-paxticle  problems 
on  the  other  hand,    although  some   success  has  been  achieved,   the  agreement  between 
theory  and  experiment  is  often  poor,  particularly  when  exchange  scattering  is   in- 
vestigated.    Moreover  the  mathematical  procedures  which  have  been  devised  to  ob- 
tain numerical  results   in  these  many-body  problems  are  not  merely  non-rigorous 

2 
but   appear  often  to  contain  internal  inconsistencies.        The   difficulties   involved 

in  deducing  the  experimentally  very  successful  Breit-Wigner  formula  for  nuclear 
reactions  are  well  known  .      Somewhat  less  well-known  are  the  difficulties  en- 
countered vrtien  attempting  to  deduce  a  variational  principle   of  the  Schwinger 
type    '-^  for  scattering  with  the  inclusion  of  exchange.     It  has  been  conjectured 
that   in  the  problem  of  the  scattering  of   electrons  by  hydrogen  atoms,   no   itera- 
tion method  stemming  from  the  usual  integral-equation  formulation  of  the   scatter- 
ing problem  can  ever  converge  to  give  exchange  scattering. 

In  our  opinion  the  difficulties  which  have  been  mentioned  stem  in  part 
from  em  incomplete  understanding  of  the  mathematical  questions  involved.      It   is 
these  questions   to  which  we  devote  our  attention.     Generally  the   solution  of  a 
scattering  problem  involves  related  difficulties   concerning  the  boundary  condi- 
tions, uniqueness,    completeness,  and  the  determination  of  a  G-reen's  function. 
There  is  a  well-developed  theory  available  for  answering  such  questions  about 
one-particle  problems  in  one  dimension.      This   theory  can  b'e  extended  to  one- 
particle  three-dimensional  problems.     The  extension  to  mfiny-particle  problems 
is  less   immediate  and,   at  this  time,   less  rigorous.     However,    the   one-particle 
problem  famishes  a  useful  and  clarifying  guide   toward  the  correct  formulation 
in  many-psu:ticle  systems. 

In  this  paper  we  shall  consider  one-particle  one-dimensional  systems. 
Many-dimensional   many-particle  cases  v;ill  be  treated  in  future  reports. 

The  time-independent   Schroedinger  equation  is 

(1)  •  ^    ^     +V(x>|r=  Bt   , 

dz 

which  we  write  in  the  form 
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(2)  -  ^  +  qi  =  X^  : 

dx 

2m7   .  _  2mE 

In' scattering  prolilems  we  seek  (for  arbitrary  positive  X)  a,   solution 
of  Eq.  (2)  which  is  continuous  and  possesses  a  continuous  derivative  over  the 
entire  interval  -00  <  x  <  00  ,  This  solution  is  to  have  the  form  of  an  inci- 
dent wave  and  a  scattered  part,  as  follows: 


l/K 


(3)  >Kx)  =  e^^'^  *+0U)  , 

where  0(x)  is  "everywhere  outgoing  at  infinity",  e.g. 

(1^)  x-»oo   :  0(i)~Ae^''^^ 

x->-oo  :  0(x)~Be"^^^  . 


We  shall  show  that  the  constants  A  and  B  are  determined  hy  the  outgoing-wave  con- 

;ic 
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dition.  The  transmission  coefficient  is  |a|   and  the  reflection  coefficient  is 


II.  Boundary  Conditions  for  a  Finite  Interval 

Bq,  (2)  is  the  special  case  p  =  1  of  the  equation 

(5)  L-|^(p^)*Ot  =  xM^ 

which  is  Itself  a  special  csise  of 

(6)  (L-X)i=0, 

where  L  is  a  linear  operator.   To  make  equation   (6)  meaningful  it  is  necessary  to 
specify  the  domain  of  L,   i.e.    the  set  of  functions  to  which  the  operation  L  may 
he  applied. 

We  shall  assume  that  the   operator  L  is  considered  in  a  finite  interval 
(x- ,x   ),   that  p(x)  is  continuous,   different iahle  and  never  zero   in  the  closed 
interval  (x- ,x   ),   and  that  q(x)   is  hounded  in  that   interval.     The  domain  of  L 
will  consist  of  those  functions  v(x)  vrtiich  have  piecewise  continuous   second 
derivatives  in  the   interval    (x^.Xp)  and  for  which  the  second  derivatives  are 
square  integrahle  over   (x^.Xg). 
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We  define  a  scalar  produce  of  fimctions  in  the  domain  as  follows: 


(7) 


[v,vfj   =     /      T    wdx  . 


Here  the  asterisk  denotes  complex  conjugate.  If  we  introdace  the  Wronsklan 
W(u-,Up)  of  two  functions  u.,u  "by  the  definition: 

(8)  WCu^.u^)  =  u^Ug  -  u^u^'  , 

then  integration  hy  parts  shows  that 


(9) 


[v.Lw]    -   [Lv.w] 


lW(v*,w) 


^ 


The  symbol 


indicates  the   difference  between  the  values  at  x  =  Zp  and 


X  =  X.    of  the  function  to  the  left   of  the  symbol. 

Equation  (9)  is  valid  even  if  q(x)  has  a  jump  discontinuity  at  a 
finite  number  of  points.      Such  discontinuities  occur,    for  example,   in  the 
quantum  mechanical  problem  of   scattering  by  a  square  well. 

Now  let  us  consider  a  subset  D  of  the   domain  of  L.     D  is  such  that 
the  fTinctions  v(x)   in  it  satisfy  the  boundary  conditions 

(10)  v(x^)  +0(v'(Xj^)  =  0  . 

(11)  vCxg)  -^-^v'U^)  =  0   , 

where  oc  and  fl   are  arbitrary  fixed  real  numbers,    includirg  possibly  *   oo    ,   and 
where  the  prime  signifies  the  derivative.     If  v  and  w  belong  to  the  subset  D, 
that   is,   if  they  both  satisfy  the  boTindary  conditions   (10)  and  (ll),   then  from 
(9)  we  find  that 

(12)  [v.LwJ   -   C^^i^I]    =  0  • 

This   expresses  the  fact  that  L  is   congjlex   self-adjoint  or  hermitlan. 

The  eigenfunctions  u(x)   of  the  operator  L  are  those  functions   in  D 
which  satisfy  equation  (5).      In  general,    for  arbitrary  values   of  X  there  vrill 
be  no  solution  of   (5)  vAiich  satisfies  both  bovindary  conditions    (lO)  and   (11 ). 
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To  find  the  eigenvalues  and  the  eigenftinctions  we  start  "by  constjructing  two 
solutions  v^   and  v     of   (5).  v-(x^,^)  satisfies  the  boundary  condition  (10), 
and  V   (x,,X)  satisfies  the  hoTindary  condition  (ll).  Then  we  shall  investigate 
for  what  values  of  X  these  two  solutions  are  essentially  the  same. 

Since  equation  (5)  is  a  differential  equation  of  the  second  order, 
a  unique  solution  can  be  specified  by  giving  the  values  of  '^  and  i|f=  at  x  =  x^ 
or  at  X  =  X  ,      To   define  a  solution  satisfying  (10),    it  is   sufficient   to  give 
\|r'(x-  )  an  arbitrary  value  and  then  to  put  iKx^ )  =  -o<\if«(x^).     However,    for 
later  purposes  it  is   important   that  the  initial  values  of  the  solution  be  in- 
dependent of  the  value  of  \.     This  requirement    can  be  satisfied  in  many  ways. 
The   following  way  is  advantageous  since  it   also  makes  sense  if  c^  =  0  or  cX  =    00: 

Put  2 

v'  (x-,X)  =     -^ —  ,  V.  (x^,\)  =  -     — 5 . 

■••        •■•  ot'^+l         ^     ^  o('^+l 

Besides  v,  (x  ,  \)   the  ftmction  Cv^{x,X),    where  C  is  independent  of  x,  is  also 

a  solution  of  (5)  which  satisfies  (10),  Conversely,  any  solution  of  (5)  which 

satisfies  (10)  must  be  of  the  form  Cv  (x,A), 

Similarly,  we  may  define  vAx,X)   uniquely  by  putting 

■/   X  ^    3+1        (       ,  >  _    B^+!S 
V '(x_,X)  =  -5 ,  V  (x  ,X)  =  - 


2  2*  2 

It  is  easily  shown  that  both  v  (x,X)  and  v  (x,X)  are  entire  functions  of  X  for 

any  fixed  vsdue  of  x. 

In  general,  v.  and  v„  are  linearly  independent  functions.  The  eigen- 
values of  L  are  those  values  of  X  for  which  v^  and  Vp  are  linearly  dependent, 
that  is, 

(13)  "^2  "  ^''l' 

because  if   (13)   is   true   thsn  Vp  is  a' solution  of  (5)  which  satisfies   (ll)  and  also 

(10).      Equation   (13 )  is  completely  equivalent  to  the  following: 

(lU)  w[;v^(x,X),V2(x,X):!   =  0 

Note  that  for  a  fixed  value  of  X,  any  two  solutions  ^^  and  ^^   of  (5) 
satisfy  the  relation 

(15)  lrCpwci'i.^2^^  =  ° 
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whether  or  not  the  boundary  conditiODa  (lO)  amd  (11 )  are  satisfied;   therefore, 
pW[^  i|/^(x,\),\y  (x,X)  ]  is  independent  of  x.  Since  p(x)  is  never  zero  in  the 
interval,  it  is  satisfied  everywhere  in  the  interval,  Conseauently,  the  zeroes 
of  W'v, ,Vp),  as  a  function  of  X,  yield  the  eigenvalues  of  the  operator 

(16)  L  =  -  fj  (p  ^)  .  , 

considered  on  the  set  D.     Note  that   because  of   (13)  an  eigenfxinction  u.  (x,\.  ) 
corresponding  to  an  eigenveuLue  X.    is  unique  except  for  a  multiplicative  factor. 

If  u(i)   is   in  D,   then  so  is  u  (x),    since  <X    and  /3     In  (10)  and  (11)  are 
both  real.      In  equation  (12)  put  v(x)  =  w(x)  =  u.(x,X   )jthen  we  find  that 

(17)  (\  -  \    )    [u^.u^H   =  0. 

* 
which  implies  that  ^^-^^i   -  0,    'that   is,  A  must  be  real.     Similarly,    if  we  put 

v(x)  =  u.(x,X.)w(x)  =  u   (x,X.  )^  vdiere  A,   f  ^i,^^  obtain  the  orthogonality  relation 

(18)  &i'^jJ    "     °* 

III,  Expansion  Theorem  for  a  Finite  Interval 

It  can  be  proved  that  the  set  of  eigenfunctions  u,  of  L  (Eq.  (l6))  is 
complete.  More  precisely,  L  has  a  denumerable  set  of  eigenvalues  X.  with  corre- 
sponding normalized  quadratically  integrable  eigenfunctions  u.  in  the  interval 
X.5  X  5  Xp  .  Also,  L  is  such  that  an  arbitrary  quadratically  integrable  function 
)|f  can  be  expauaded  in  terms  of  u.  as  follows: 

(19)  ^  =  11 


a.u 


l"i* 

Here 

(20)  Cu^.Uj:  =  6^j 

and  then 

(21)  a^  =  Cu^.t:]  . 

It  must  be  noted  that  Eq.  (19 )  is  true  only  in  the  sense  that 

(22)  liiD   U  -  E  Vn'  ^  -  ^  Vj  '   °- 
n— >  00  1_     n  J 

For  more  restricted  classes  of  functions   the  series    (19)  may  converge  everywhere, 
but  not  necessarily  to  the  value  of  ijf,  \*iile   for  more  restricted  classes   still, 
it  may  converge  to   the  value   of  i|f  at  all  points   in  the   Interval  Xt    f  x  £  Xp. 
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With  the  help  of  the  following  Schwarz  inequality  for  any  two  quadratic ally 
integratle  functions  f,g: 

(23)  [f.g]^     <  [f.f]    [g.g]. 

we  find  that 

i2h)  lim     Cf.   Ylwl   =   &.^D    . 

n.-»oo  n 

since 

1  1 

(25)  lim     Hf.i:  a^u^-M^:]    <  lim     [f  .f]^  [i^  -  ^  a^u^.  i|r  -  ^  a^u   1  ^     . 
n-»ao         n  n->oo  >-  n  n  -' 

Thus  for  purposes  of  integration,   the  differences  "between  ^  and  JZ  ^  ''^     ^'^y  ^® 
ignored. 

From  (19)  azid  (21)  we  conclude  that  we  may  write 

(26)  H^A^^  ii4*(3c«(  =  6(x-x«) 

i 

since  for  any  i 

^2 

(27)  lim     /djc'  5~u  (x)u*(x«)  ^Kx')  =  lim     5"  a  u  (x)  =  /   djc6(x-x»  )iKx' )   . 

J  *—    n         n  V        —    n  n  -^ 

n->oo  n  n-»oo   n  x^ 

It  can  also  he  demonstrated  directly  from  the  properties   of  the  Green's   function 

that 

*. 


(28)  lim     y-     rdx^Kx}u  (x)  u^(x')  =  \|f(x')  . 

The  Green's  function  G(x,y.',\)  is  a  solution  of  the  equation 


(29)  (   -I^PST'*"'!"^^  G(x,x«.>J  =  6(x-x') 

which  satisfies  the  "boundary  conditions   (7)  and  (S).     Eq, (29)  Implies  that  G(x,x',X) 

has  a  discontinuity  at  x  =  x'   such  that 

(30)  [fe=u.^'.x)]   '      =-  ^  . 

x«-0 
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Eemem'bering  (11)  we  see  that  for  X  f  \.  ,   G(x,x',X)  can  be  written  as 


X  <  x'  :  G(x,x«,\)  =  - 

(31) 

X  >  x«  :  G(x,x«,\)  =  - 


vv{y^,r^) 


pW(v^,V2) 


The  functions  v^  and  Vp  are  defined  as  in  section  II. 

For  X   =  X.,   en  eigenvalue  of  L,  the  zero  of  pW(v  v„)  invalidates 
Eq.  (3l)»   This  is  to  be  expected  since  the  homogeneous  equation  correspond- 
ing to  Eq.  (29)  has  a  non-zero  solution  and  therefore  i\  =  A.  .   On  the  other 
hand  if  X  ?^  X.,      the  Green's  function  exists  and  is  unique.   Its  existence 
when  X  ^  X.    is  insured  by  the  explicit  construction  of  G  in  Eq.  (31),  and 
the  uniqueness  of  this  Green's  function  follows  because  if  G^  and  G  aire  two 
solutions  of  Eq.  (29)  satisfying  the  boundary  conditions,  then  G  =  6  -G 
satisfies  the  boundary  conditions  and  the  equation 

(32)  ^  -  lrP|^*  q-X)G(x,x',X)  =  0. 

We  have  already  shown  that  for  X  ^  X.    there  are  nc  non-zero  solutions   to  Ea. 
(32)   (which  is   equivalent  to   Eq.    (5))  satisfying  the  boundary  conditions. 
Hence  the  Green's  function  is  unique. 

Combining  Eq.    (26)  with  Eq,    (29)  we  obtain  the  following    well-kr-own 
expansion  of  the  Green's  function,  valid  for  X  =  X.i 

u^(x)  u*(x') 


(33)  G(x.x',>v)  =     ~       ^  ^ 


Eq.  (31)  shows  that  G(x,x',X)  =  G(x',x,X).   This  symmetry  is  not  manifested  in 
2q.  (33),  ^t  it  is  implied,  since  in  the  problem  we  have  been  considering, 
u.(x)  differs  from  a  real  fimction  by  a  complex  multiplicative  factor  only, 
and  u  (x)u*(x')  is  real.  For  the  same  reason  the  asterisk  may  be  attached  at 
will  to  either  u(x)  or  u(x')  in  Equations  (26)  through  (28). 

We  have  thus  far  taken  X   to  be  real.   However,  we  can  consider  solu- 
tions of  Eq.  (5)  for  complex  X,   and  as  we  have  already  pointed  out  it  can  be 
shown  that  the  functions  y-^{x,X)   and  v„(i,X)  which  satisfy  the  boundary  condi- 
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tlons  (7)  and  (8)  respectively  and  which  are  normalized  independently  are  analytic 
functions  of  X.     Consequently  G(x,x',\)  as  defined  hy  Eqs.  (29)-(5l)  is  an  ana- 
lytic function  of  A  everywhere  ill  the  complex  X-plane  except  at  eigenvalues 
X  =  X   where  G(x,x',X)  has  a  pole.  The  right  side  of  Eq.  (33)  is  similarly 
analytic,  and  hence  we  see  ty  analytic  continuation  that  Eq.  (33)  remains 
valid  for  complex  X,  Moreover 

^'^^^  ~  2ni  /*  <3AG(x,y',X)  =  XI  u^(x)uj(x')  =  6(x-x'). 

C  i 

In  Eq.  (3^),  C  is  a  positive  (counter-clocl:wise )  contour  enclosing  all  the  poles 
of  G(x,x',X). 

IV,  Expansion  Theorem  for  More  General  Operators 

Equation  (3^)  shows  that  the  integration  of  the  Green's  function  over 

a  suitable  contour  in  the  complex  X-plane  yields  the  complete  set  of  normalized 

eigenfunctions  of  L.   This  result  cnn  "be  applied  to  more  general  types  of  opera- 

g 
tors,  for  example  those  having  continuous  spectra.  Titchmarsh  has  given  a 

rigorous  treatment  of  such  operators  using  the  theory  of  the  Green's  f\xnction, 

A  summary  of  the  theory,  with  applications  to  electromagnetic  problems,  has 

9 
teen  given  by  Marcuvitz-^. 

In  what  follows  we  shall  discuss  the  method  and  give  a  few  elementary 

examples  by  way  of  illustration.  For  a  complete  mathematical  justification  the 

g 
reader  is  referred  to  Titchmarsh  . 

We  start  with  the  differential  equation  defining  the  Green's  function, 


namely 

(L  -  X)G(x,x',X)  =  6(x-x')  , 
where  L  is  some  linear  differential  operator  defined  over  some  domain,  and  X  is 
an  arbitrary  complex  number.  The  domain  of  L  is  to  be  chosen  that  L  is  complex 
self-adjoint. 

Consider  the  integral  of  G(x,x',X)  over  a  circle  of  radius  E  in  the 
complex  X-plane.   We  put 

(35)  Ir=  2Fi  /rGC^.^',^)  ^' 
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Eq.  (3^)  shows  that  for  the  relatively  simple  operators  considered  so  far, 
whose  eigenfunctions  form  a  complete  set, 

(36)  lim   Ip  =  -6(x-x')  . 

R-»oo 

For  more  general   operators  we  must   investigate   the  "behavior  of  the  Green's 
function  for  large  values  of  A  and  prove   that    (36)   is   correct.      The  Green's 
function  continues   to  he   constructed  from  those  solutions  of 

(37)  (L  -  X)  u  =  0 

which  satisfy  the  boundary  condition  at  one  or  the  other  end  of  the  interval. 
By  examining  the  behavior  of  these  solutions  for  large  A,  the  validity  of  Bq. 
(36)  can  be  established.   For  example,  if  we  put  p  =  1  in  Eq.  (I6)  then  the 

solutions  u(x,A)  generally  behave  like  e"  for  large  |a| ,  and  we  have  the 
approximate  relation 

i/A  |x-x'| 


(32)  G(x,x')  =  - 


2 
If  we  now  prut  A  =  k  ,  we  obtain 


e 


21/A 


lim  Ij^  =  lim  -  ^^     -— — dA  =  -  -  y  dke   I     I  =  -6(x-x')  . 


Consequently,  londer  these  circ\imstances  Eq.  (36)  has  been  established,  and  it 
is  evidently  quite  generally  true  that  Eq.  (36)  is  valid  if  Eq.  (38)  holds  for 
large  |a|  . 

Once  (36)  has  been  established,  the  circle  of  radius  R  csm  be  modified 
into  any  suitable  contour  C.   If  G(x,x',A)  has  only  poles,  then  G  can  be  taken 
as  a  set  of  small  circles,  one  around  each  pole,  amd  thus  by  evaluating  the  resi- 
dues we  obtain  the  expansion  of  6(x-x')  as  a  sum  of  discrete  eigenfunctions.  If 
G(x,x',A)  also  has  branch-point  singularities,  then  C  can  be  taken  as  the  set  of 
small  circles  around  each  pole  plu?  contours  around  each  of  the  branch  cuts.   From 
this  contour  we  obtain  the  expansion  of  5(x-x')  as  a  sum  of  discrete  eigenfunc- 
tions (residues  from  the  poles)  plus  an  integral  of  continuum  eigenfunctions  (the 
integrals  over  the  cuts).  These  details  will  become  clearer  when  the  illustrations 
of  the  next  section  are  considered. 
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From  the  expansion  of  the  6  ftinction,  the  expansion  of  any  arMtrary 
function  csm  he  obtained  as  follows.  Suppose  that 

6(x-x')  =  YZ    ^i^^^  \^^*h 
then,  multiplying  hy  f(x')  and  integrating,  we  ohtaln 


f (x)  =  Yl  u^(x)yf(x»)  u*(x')dx'. 


V  .  5xf:)an8ion  Theorem.   Illustrations  for  a  Finite  Interval 

We  present  some  examples  in  which  the  theory  of  the  preceding  section 

Is  us'^d. 

Choose  p  =  1,  q  =  0  in  Eq.  (5);  then 

.2 
(39)  L  =  -  ^  . 

dx 

Suppose  that  the  interval  is  0  <  x  ^  ti.  The  spectrum  will  of  course  depend  on  the 
boundary  conditions.  Thus  if  v(o)  =  ri'n)   =  0  are  the  boundary  conditions,  we  have, 
for  arbitrary  complex  \ 

(liO)  y^iT,\)   =  9in\^  X 

CjI)  Vg^^'^^  ""  sinVx  (x-tt) 

and,  substituting  in  T^q.  ('^l), 

X  <  X'  :  G(x.x'.X)  =  -  «inxVXsin  (x '-tt) /T 

•/x  sin  TT  yX 
(^■2) 

x>x«  :  c;(x,x'.x)  =  -   ^^"  (^-WV^«i"x'/r    ^ 

■/x  sin  TT  -/T 
Note  that  for  large  values  of  X,0  <  arg  X  <  2n,G(x,x',X)  behaves  like 

^i/X  Ix-x'l  /2i/X  ;   consequently  (36)  can  be  proved  valid.  The  poles  of 

o 
G(x,x',X)  are  at  X  =  n  ,  n  being  any  positive  integer.  There  is  no  pole  at 

X  =  0,  where  G  is  finite,  nor  is  there  a  branch  point  at  X  =  0,  since  although 
/X  occurs  in  the  numerator  and  denominator  of  G,  G(x,x»,X)  returns  to  its  origi- 
nal value  as  X  circles  the  origin.  The  contour  C  can  be  deformed  to  become  the 

contour  C  (Flg.l),  so  for  arbitrary  x,x'  the  evaluation  of  the  integral  (35) 

2 
simply  involves  calculating  the  residues  of  G  at  the  poles  X'  =  n  . 
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c 


X^v_ 


c 


Fig.   1 
Contour  for  G(x,y.',\)   of  Ecj.    (i)2). 

Using  Cj(x,x',\)  for  x  <  x\   we  o"btain  the  following: 

/i,-\            1      f  j-,r  /        l^^        V2  sin  nxsin  n(x'-TT)       v—     2      .  .  , 

(^3)     -  -rn- /  <i>^&(x.x',X)  =  2_  Z  ^^„  „„ =2—    -  sin  nx  sin  nx'. 

2TT1  /  „  '—  TT  cos    nTT  — ^,     TT 

^  C  n  n=l 

The  normalized  eigenftmctions  are  therefore 

(^^)  n^M   =yf^8in  nx, 

and  in  the  interval    (O.rr)  an  arTjitrary  square   Integrable  function  f(x)  hfis  the  expansion 


(M) 


2  JO  / 

f (x)  =  —  )       sin  nx  /  dx'   f(x')   sin  nx', 

n=i  0 


This  is  the  well-lcnown  expansion  in  a  series  of  sine  functions. 

If  we  choose  the  boundary  conditions  ▼•(o)  =  ▼•(tt)  =  0,  then 


(^6) 


V,  (x,>  )  =  cos  yX  X 
▼pCx.X)  =  cos  v/X    (x-tt) 
W(v^,v,)  =  A  sin  TT-Zx 


X  <  x»    :    G(x,x',X)  =  - 


cosx  yX  co8(x'-tt)  Vx 
•/x  sin  -^  /\ 

X  >  X'    :    G(x.x',X)  =  -  C03(x-tt)  VX  cos   x '/x 

V^     sin   TT  y/x" 
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The  poles  are  at  X  =  n  and  also  at  \  =  0.  G(x,x',x)  is  single  valued, 
and  the  conto^or  is  the  same  as  in  Fig.  1  above,  except  that  the  residue  at  \  =  0 
must  also  'be   evaluated.  Using  Eq,  {h6)   we  have 

{t'?)  -   -i-  /  dX  G(x,x',X)  =  XI   -  cos  nx  cos  nx'  +  - 

2"i  '/C  n=l 

which  yields  the  followlnfr  cosine  expansion  in  the  interval  (O.fr): 

,  TT  00.        p 

(L8)  f(x)  =  -/  dx'f(x')  +  ^  H  cos  nx  /  dx'f(x')  cos  nx'. 

"''  0  ^  n=l       "   0 

Eqs.  (^"7)  and  (^8)  can  of  course  he  comhined  in  the  familiar  way  to 
ohtaln  the  expansion  of  an  arhitrary  function  in  the  interval  (-tt,tt)  in  a  series 
of  sines  and  cosines.  To  ohtain  this  expansion  directly,  it  is  necessary  to  im- 
pose periodic  "boundnry  conditions: 

u(-tt)  =  u(tt) 
^^"^  u'(-tt)  =  u'(tt). 

These  so-called  mixed  "boundary  conditions  are  not  of  the  form  (?)  and  (8).  How- 
ever, it  is  easy  to  see  from  Eq.  (lO)  that  Eq.  (i+9)  still  guarantees  Eq.  (12). 

d2 

Consequently  the  operator  L  =  -  — r-  remains  hermitian  self-adjoint. 

dx^ 

We  can  find  the  expansion  theorem  for  L  in  the  domain  defined  hy  (^9) 
by  the  method  of  the  preceding  section.  Since  independent  soltitions  of  the  differ- 
ential equation 

u."+  Xu  =  0 

are  sin  y'X  x  and  cos  ■/x  x,  we  may  assume  that 

(50)      G(x,x',X)  =  a  sin  ^  x  +  h  cos  -/x  y,  x  <  x' 

=  c  sin  yX  x  +  d  cos  yX  y,  x  >  x'  . 

The  coefficients  a,  h,  c,  d,  are  determined  hy  the  boundary  conditions  {k^)   and 
by  equntlon  (?0).  We  find  that 

nt        I  ^  N     cos  vX  X  cosVxCtt-x')  -  sinitX  X  sinvX(TT-y')        , 
G(x,x',X;  =  -  * ' — ^ ' " ' — ^ '-   ,  X  <  x' 

2  /x  sin  TT  y/x" 

(■51)  r-  r-  r-  r- 

cos  yX  X  cos  VX(th-x'  )  +  sin  yX  x  sin  VX(th-i«)        , 

—    —  ,    X   ^  X    • 

2  yfk  Sin  tt  yX 

The  contour  and  the  singularities  of  the  Green's  function  are  the  same  as   in  the 
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case  of  the  cosine  series  above.   We  obtain 

00 


1   /"  ^ 

-   T— 7  /d\G(x,x',X)  =  ~  X_  (cos  nx   cos  nx'  +  sin  nx   sin 


nx 


)*i: 


n=l  2" 


and  this  Implies  that 


=  ^  /  dx'f  (x')+  ^XI  |cos  ni  /  dx'f  (x')co8  nx'  +  sin  ni  /  dx«f(x«)sln  nx«  [  . 

This  example  1%  very  helpful  In  illuBtratlng  what  occurs  with  the  Green's 
function  in  more  than  one  dimeneion.  We  see  that  G(x,x',\)  may  not  be  of  the 
simple  form  ^iven  by  (3l).  namely  a  function  of  x  multiplied  by  a  function  of  x' 
and  divided  by  a  function  of  X..   Instead,  as  le  shown  in  Eq.  (51),  the  Green's 
function  may  be  a  sum  of  products  of  functions  of  x  and  functions  of  i',  all  di- 
vided by  a  function  of  X.   In  such  a  case  each  term  in  the  sum  will  become  an 

eigenfunction  and  we  will  have  degeneracy.   In  the  present  case,  to  each  eigen- 

2 
value  X  =  n  there  correspond  two  eigenfunctions  sin  nx  and  cos  nx. 

VI, Boundary  Conditions  for  an  Infinite  Interval 

When  p(x)  becomes  zero  or  q(x)  becomes  infinite  at  the  endpoints,  or 
when  the  interval  le  Infinite,  the  endpoints  are  termed  singular,  and  the  bounfl?^ry 
conditions  must  be  examined  more  closely.  The  problems  of  interest  in  quantum 
mechanics  ar©  almost  always  singular  in  this  sense.  The  necessity  for  reformu- 
lation is  apparent  from  the  fact  that  in  the  singular  case,  if  u(x)  is  a  solution 
of  the  differential  equation  (5),  the  limits  of  u(x)  or  of  Ti'(x)  generally  do  not 
exist  at  the  endpoints  of  the  interval.   Consequently,  the  boundary  conditions 
('')  are  not  applicable.   In  partic^ilar,  it  is  known  that  sometimes  the  mere  re- 
quirement of  quadratic  integrabllity  is  a  sufficient  bcandary  condition.  This 

10 

was  elucidated  by  Weyl 

It  turns  out  that  the  bcmdary  conditions  at  singular  endpoints  are  best 
understood  by  considering  the  eigenvalue  problem  for  complex  values  of  X,  as  we 
did  in  the  previous  sections.  The  essential  results  are  as  follows'.  Suppose 
p(t)  is  continuous  and  greater  than  zero  everywhere  in  the  open  interval  x,  <  x  <  x_. 
If  x^  =  -CD  or  Xp  =00,  or  if  p(x)  vanishes  at  x  =  x,  or  nt  x  =  x-,  then  x  =  x^ 
and  X  =  x»  are  singular  points  of  the  differential  equation. 

Let  X-  be  any  point  in  the  interval  (x-,Xp)  and  let  w,(x,X)  and  w-(x,X) 
be  Solutions  of  Eq.  (5)  satisfying  the  conditions  w^(x_)  =  1,  w'(x_)  =  0  and 
W2(x-)  =  0,  w^(x^)  =  1  respectively. 
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Consider  those  solutions  of  Eq.  (5)  which  satisfy  the  condition 
u'(x'.X)  +  Y  ^(^'.>^)  =  0. 
where  x'  is  some  fixed  point  between  x-,  and  x  .  Since  the  general  solution  of 
Equation  (5),  except  for  a  constant  factor,  is 

u(x,X)  =  v^{x,X)   +  »?  w^Cx.X)  , 

we  have 

w_[(x«.A)  +  rjWgCx'.X)  +YWj^(x',X)  +  rr^v^{j:t^x)   =  0  . 

From  this  equation  it  follows  that 

rw  (x',\)  +  w'(x',X) 

(52)  n  =  — 7 • 

rw^Cx'.X)  +  WgCx'.X) 

If  T  goes  over  all  real  values  from  -00  to  +00,  then  the  values  of  rj^  given  "by 
Eq.  (52)  will,  for  fixed  x',  describe  a  circle  in  the  complex  «^ -plane, 

Weyl  showed  that  if  X   is  not  real, the  circle  defined  ty  Eq.  (52)  for 
x'  =  X.  will  contain  the  circle  defined  for  x'  =  x  if  x^^  >  x  .  As  x'  approaches 
x^  the  circles  either  shrink  down  to  a  limit  circle  or  to  a  limit  point.   In  the 

first  case  it  can  "be  shown  that  for  every  real  or  complex  value  of  «  ,  we  have 


/^2  2 

dx  Ju(x,X)j  <  00  . 


^3 
In  the  second  case  there  will  exist  one  and  only  one  value  of  i^  such  that  Eq, 
(53)  is  satisfied. 

Since  at  an  endpoint  in  the  limit-point  case  there  is  only  one  quad- 
ratically  integrahle  solution  of  the  differential  equation,  we  need  not  and  in 
fact  can  not  impose  a  'boundary  condition  on  the  solution  there.  The  mere  re- 
quirement of  qiiadratic  integrahility  determines  the  solution  iiniquely.  Note 
that  these  statements  are  valid  only  for  complex  values  of  X;  when  X  is  real 
there  may  "be  no  solutions  which  are  quadratically  integrahle.   For  example, 
consider  the  equation 

(51+)  u"  +  Xu=  0 
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In  the  interval    (-00,    00).      If  the  imaginary  part  of  X  is  positive,    then  e^'^'^^  is 

quadratically   Integrahle  at  +00,   while  e"^*    ^     is   quadratically  Integrable  at    -oj. 
Since  no  other  solutions  are  quadratically  integrahle,    it   is    clear  that  ^00  are 
limit-point   cases.      When  X  is  real  there   is   no  q\iadi-atically   integrahle   solution 
of   (5^+). 

On  the  other  hand,    in  the  limit-circle  case  there  are  an  infinite  num- 
"ber  of  solutions  which  are  quadratically   integrable.      In  order  to  specify  the 
domain  on  which  L  acts  so  as  to  make  L  hermitian,    it  is  absolutely  necessary  to 
specify  an  appropriate  boundary  condition  at   the  endpoints.      Once  this   is   done, 
the  well-known  theorems  about  the  spectrum  of  the  operator  can  be  used,   sind  the 
orthogonality  of  the  eigenfunctions  follows  automatically.      These  considerations 
would  have  greatly  simplified  some  recent  work  by  Case       on  singular  potentials. 

Ve  may  summarize  the  discussion  of  this  section  as   follows:      In  the 
llBit-point  case,    the  mere  requirement  of  quadratic  inteecrability  will  determine 
uniquely  the  solutions  and  the  Green's  function  for  complex  values  of  X.     By 
analytic  continuation  we  will  get  the  solutions  and  the  Green's  function  for 
real  values  of  X.      In  the  limit-circle  case  boundary  conditions  must  be  given 
at  both  endpoints,   and  then  we  proceed  as   in  the  case  of  the   finite  interval, 

VII.   Expansion  Theorem.      Illustrations  for  an  Infinite   Interval 

We  now  illustrate  the  manner  in  which  the  considerations  of  Section  V 

are  combined  with  Eq.    (36)  to   obtain  the  normalized  eigenfunctions  when  the 

spectrum  is  continuous. 

A.    Fourier  integral. 

.2 

Consider  again  the  operator  L  = s-  ,   over  the   interval    -00    <   x  <    00. 

dx'^ 

Both  endpoints  are  now  singular.  A  general  solution  of  Eq.  (5)  for  arbitrary  com- 
plex X  is 

(55)  u(x,X)  =  e^^''  -e   re"^'^^ 

(see  Eq,    (52)).      Define  ^  so  that   it   is  real  and  positive  on  the  positive  real 

X  axis.     This    implies  that   if  X  =    |x|e^®   ,    then  X^/^  =    |x|-^^^  e^®/^,   0  i  0  <  2n. 
For  non-real  X,    Im  X  >  0,  and  the  only  quadratically  integrable  solution  as 
x->oo  is  T  =  0  in   (55)  above.      In  other  words,   we  have  a  limit  point  at  x  =    00, 
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There  is  similarly  a  limit  point  at  x  =  -oo,  the  only  permitted  solution  at  that 
endpoint  being  u(x,A)  =  e~  '  ^  .  We  have  therefore 


Ut^(x,X)  =  e~  ^^  ^  , 
u  (x,X)  =  e  *^  ^  , 


(56)  -2 

WCu^.u^)  =  21/X  , 
and  then  the  Green's  function  is 


(57) 


X  <  x«    :    G(x,x',X)  =  - 
X  >  x»    :    Ou,x»,X)  =  - 


e-i  A: 

ic     i 

e 

/K  X* 

2i 
^i/Xx 

75^ 

-i 

e 

7X  x' 

21  yx" 


The  Wronskian  has  no  poles,  hut  there  is  a  hranch-point   singularity 
at  X  =  0     and  consecjiently  G(x,x',\)   is  not  a  single-valued  fvmction  of  X. 
Since  we  have  defined  X  in  the   range  0  <  arg  X   <  2n  only,  we  may  draw  a  cut 
eilong  the  positive  real  X  axis  and  then  we  have,  when  x  <  x'. 


(5S) 


-  ^  f  ^  ^i-.-' ^)  -    ~I 


dA 


iJX  (x'-x) 


"C  -'^    "C  21^" 

with   the  contour  as   in  Fig.    2.      Explicitly  calculating  the   integrals  along 


////////////////////// 


Figure  2 
Contour  for  Eq.  (58) 


Cross-hatching  shows  cut  along  positive  real  axis. 
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the  upper  and  lower  aides  of  the  cut,  we  obtain 

/   X     1   /*  1    9P    „i  v^  (x'-x)   ,   00     -i  A"  (x«-x) 

(59)  -5Sr/««-.-.M--4/dA  2 ,1^^  dx'  '^'"    . 

On  the  right  side  of  Eq.  (59),  X  ie  real  and  positire.   We  infer  that  the  spectnun 
is  continuous  and  nms  from  0  to  00,  and  that  there  are  two  eigenfunctions  for 
each  value  of  X,  namely 

(60) 

The  functions  u^  and  u^  of  Eq.  (60)  are  properly  normalized.   The  cor- 
responding expansion  theorem  is 


-CO   i  y\  X   .00  4  /T  I 

f(x)  =  ^/dx2_^  /   dx.  f(x.)e-^>^^' 

(61) 

,   .00   -i  yX  X  -00  .  rr^f 

^   i_y  dX  ^-— /   dx'  f(x')e^  y^^   , 


If  we  introduce  k  =  /x    ,  Eq.    (61 )   takes   the  more  familisr  form 


°°  --     '°°  -ikx' 

e 


(62)  f(x)  =  ^/      dk  e*"^  /      dx'  f  (x') 

'       -  CO  -  00 

Because  of  the  depeneracy,  different  sets  of  eigenfunctions  exist 
which  leed  to  nlternptive  forms  of  the  expansion  theorem.  Thus  Eq.  (59)  could 
"be  written  rs 

1   /"a^    r(        I  ^^    1  /^i  cos  yr  (x'-x) 
-^  Jc^''   G(x.x'.X)  =  ^y  dX  — ^^ 

(63) 

=  — -  /       —   Fcosx  /x  cosx'  i/\  +   sin:>  Jx  sinx'  Jxj 
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Consequently  another  set  of  normalized  eigenfanctions  is 


/  , N    1   sinT  y\ 
u  U,X)  =  -—  — -fjf- 


(6^) 

y^  (X) 


f     >  \    1  cosx  yX 


and  the  corresponding  expension  theorem  is 

0    /x     "^00 
(65) 

.00   .,.__^   .00 


^l.f^^22S^jLL.f      dx'  f{x')co8x'^  . 


If  we  put  X  ■  k^^  then 


dk  / 

0       -  00 

Suppose  the  interval  is  0  <  y  <oo.   There  is  still  a  limit  point  nt 


(66)  fix)    =  ^  /   dkA    djc'  f(x')  C08Tr(x-T'). 


+  00,  "but  at  X  =  0  we  have  a  limit-circle  cnse.   We  impose  the  honndnry  condition 
u=Oetx=0.  We  have  now 


Ti- (t,X)  =  sin  /k  X 

f     .\         i  /x  X 
u  (x,X)  =  e 

V{M^,Vi^)    =  -  /X 


^^"^^  X  <  X'  :  G(x.x'.X)  =  sinx  /x  e^  ^^ 


X  >x'  :  C}(x,x',X)  =  - 


^>^^sinx'/x 


(68) 


.  a_y ,,,(,.,.,,)..  a_M  sin  x^e^/^^'. 


c/^ 
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The  contour,   the  cut,   and  the  definition  of  /X  are  as  in  Fig.    2  and 
Eq.    (58).      From  Eq.    (62)  after   integrating  along  the  upper  and  the  lower  sides 
of  the  cut,  we  have 

^  00      VA  0    /X 

=  ^  y  —  sin  X  yx  sin  x'  yx" . 

o   yx 

In  the  last  equality  of  (69),  X  is  real  and  positive.   Eq.  (69)  shows  that  the 
spectrum  extends  from  0  to  00,  and  that  there  is  only  one  eigenfunction  for  each 
A,  namely 

(70)  u(x,X)  =  -i-  °^^,\^   . 

Using  the  expansion  theorem,  we  find  that  in  the  interval  0  £  x  <  oo 

(71)  nx)   =  i-  /°°dA  "^^^^   /°°dx'  f(x')  sin  X'  /K 


2 
or,  if  we  put  X  =  k  , 


0       ^ 


00  00 


(72)  f(x)  =  |y  dk  sin  kxy  dx«  f(x»)  sin  kx*  . 

0  o 

B.  Alternative  representations  of  outgoing  Green's  function. 

In  order  to  "bring  out  some  important  points  we  shall  ohtain  some 
alternative  representations  for  the  Green's  function.  Consider  the  following 
integral  over  a  circle  of  radius  R: 

^'^'  R    2ni  J^        X'-X   • 

where  G(x,x',X')   is  the  Green's  function.      V/e  assume  that  as  R  approaches   in- 
finity,  G„  approaches  zero;      then  evaluating  G^  hy  residues     we  have 

where  C   is  a  contour  which  encloses  all   the   singularities,    poles,    or  "branch  cuts 
of  the  Green's  function  hut  does  not  enclose   the  point  X.      As  an  illustration. 
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consider  the  Green's  function  given  ty  (57).  We  have 

i/A'  |x'-x| 


(75) 


G(x.x',X)  =  -T^/  dX-  ^ 


c    (^'-^)/^• 


where  C  is  a  contour  enclosing  the  "branch  cut.   In  order  that  G^  approach  zero 
as  R  approaches  infinity,  the  cut  must  "be  chosen  so  that  the  imaginary  part  of 
/\'  is  positive.   This  implies  that  the  cut  must  he  taken  along  the  positive 
real  axis  and  that  0  <  arg  \'  <  2n  .  In  other  words,  the  contour  C  in  the  com- 
T)lex  \'-plane  is  precisely  the  contour  of  Fig.  2.  Put  X'  =  k  in  (75)  and  it 


becomes 

QQ    I-  ..  I   .   1         ..  I  _t  _i  -I        00 

cosk  (x'-x) 


°°   •"iklx'-xl      -iklx'-xP      °° 


^       fT  ^lk|X'-X|        g-lK|X'-X|       ^        r 

(76)  G(x.x-,A)=  2^y  dk  -^ +  ~T2-—  =2?ry. 


k  -\  k*"-  A 


dk 


-'-oo    (k-yA)(k+yX') 


Equation  (76)  is  valid  for  all  X  such  that  0  <  arg  A  <  2n.  Consequently, 
the  pole  at  k  =  /X"  is  in  the  upper  half  of  the  complex  k-plane  while  the  pole  at 
k  =  -/X  is  in  the  lower  half  of  that  plane.  Note  that  k  runs  along  the  real 
axis  by  virtue  of  the  fact  that  we  pat  A'  =  k  in  (75),  where  A'  ran  along  the 
real  positive  axis  on  either  side  of  the  cut.  However,  the  contour  in  the  k- 
plane  may  he  deformed  provided  no  singularities  are  crossed.  By  deforming  the 
contour  dovmward  to  the  ri^t  of  the  origin,  and  upward  to  the  left  of  the  origin, 
it  is  seen  that  the  last  integral  in  Eq.  (76)  defines  an  analytic  function  of  A 
as  the  iaiaginary  part  of  A  approaches  zero  for  positive  A.   Consequently,  for 
real  A  >  0,  the  contour  in  the  k-plane  must  be  the  one  shown  in  Fig.  3,  this 
being  the  only  one  for  which  Eq.  (76)  is  valid  when  G(x,x',A)  is  the  Green's 
function  of  Eq.  (57).  This  Green's  function  would  ordinarily  be  termed  the  out- 
going Green's  function,  so  we  have  obtained  an  integral  representation,  with 
correctly  chosen  contour,  of  the  outgoing  Green's  function  for  the  operator 

dx 


-I     •     *  • 1  •  r 

0  \^ 


->rr 


Fig.  3 
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There  Is  another  interesting  point  ahont  the  representation  of 
G(t,x',X)  given   in   ("^6).      The   Isat   intef^ral  contains  cos  kCx-rx')  and   this   seems 

to   indicate  that  G(x,t',\)   is  mode  up  of  ontgoiT\e  waves  e  and  also  in- 

-lk(x-x') 
coming  waves  e  .     However,  we  can  easily  see   that    the  incoming;  part  is 


Tiero  when  we  evaluate 


1.     r 


-ik(x-x') 
^  dk. 


-  *       k^-X 


For  X  approaching  oo,  the  contoiir  may  he  closed  in  the  lower  half -plane,  and  then 
the  integral  hecoraes      .  r-  ,        , v 

an  outgoing  term.  This  result  demonstrates  that  by  merely  glancing  at  an  inte- 
gral it  is  not  possihle  to  determine  its  asymptotic  behavior,  or  whether  it  is 
made  up  of  outgoing  or  incoming  waves.   This  statement  is  obvious,  but  its  impli. 
cations  have  not  always  been  understood,  especially  in  many-|Brticle  problems 
where  different  representations  of  the  Green's  function  are  possible. 

VTTI.  Transmission  and  Heflection  Coefficients  in  One  Dimension 

We  ret\irn  now  to  the  discussion  of  scattering  in  one  dimension.  We  are 
Interested  in  a  solution  of  the  equation 

(77)  _  1^  +  qw  =  X  „, 

dx'^ 

for  r°al  X.  The  solution  \|/(x)  should  be  of  the  form 


(78)  vv(x)  =  e^  ""^  ^+  ?5(x). 


-  i/X 


where  0(x)  is  everywhere  outgoing  at  infinity,  that  is, 

(79) 


0(x)  'v  A  e*^""  ,  X  -^00 


0(x)  -^B  e"^  ^^  ,  T  ->-oo. 

The  physical  interpretation  of  the  problem  Is  this:  a  plane  wave  e 
cones  In  from  -oo,  is  scattered  by  the  potential  q(x),  end  thereby  causes  the 
appearance  of  an  outgoing  wave  A  e  "^    at  +oo,  and  an  outgoing  wave  3  e~  '^ 
at  -00.   It  should  be  realized  that  the  assiojnptiona  (78)  and  (79)  about  the  be- 
havior of  'i'(x)  are  valid  only  if  q(x)  is  not  too  large  at  infinity.  For  example, 
if  q(x)  were  the  Coiilomb  potential  l/x,  then  (''8)  and  (''9)  would  be  Impossible, 
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Sufficient  conditions  for  the  validity  of  ('^8)  and  (79)  are  th^t  q(x)  vanish  out- 
side some  finite  interval  or,  more  generally,  that 

.00 


(80)  J     lq(x) 


dx  <oo, 

-00 


We  assume  (80 )  Pnd  consider  what  happens  when  a  plane  wave  e       comes 
in  from  minus  infinity.  Some  of  it  will  he  transmitted  to  plus  infinity,  while 
the  rest  of  it  will  "be  reflected  hack  to  minus  Infinity.   Let  u^(x)  represent 
the  solution  of  (77)  which  has  these  properties.  Since  we  are  interested  only 
in  the  stationnry  case,  we  may  write  the  asymptotic  relations 

U2(x)  -  e^  ^^  -  R_  e-^*^^  .  x  -^-oo 

(81)  u^(0-T_e^>^^         .x-$^co. 

where  E  is  the  reflection  amplitude  coefficient  and- T  is  the  transmission 
amplitude  coefficient.  Because  (77)  is  a  read  equation,  the  conjugate  to 
Up(x),  namely  u„(x),  is  also  a  solution  of  (77),  ai"i  from  (81)  we  have 

u(x;-~e  '     +R  e  *   ,  x  -^ -oo 
(82) 

u*(x)  -^  T*  e"^  '^  '         ,  d  ^  00  . 

We  have  mentioned  before  that   the  Wronsklan  of  any  two  solutions  of 
(77)  is  constant   for  all  values  of  x.     Let  us  apply  this  result  to  the  solu- 
tions U2(x)  and  Up(x),     We  have 


w|:u*(x),U2(x)3        =  U2(x)u2(x)  -  u^ {x)^x^{li)  «  +21  A  (1-|e_I^) 

-  00 

-00 

wru*(x).u_(x)i  =  +21/5:  |t_i^  . 

'^  '^  +00 

and  therefore 

(83)  1-   lHj2=   1tJ2  ,  or   |Rj2+  |Tj2=  i. 

Thus  the  sum  of  the  reflection  and  transmission  is  equal  to  unity.  Note  that  the 
quantum-mechanical  current  is  defined  as  ^-7-  WvUp.UpJ,  and  that  because  of  the 

properties  of  the  Wronsklan  this  current  is  independent  of  x. 
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Similar  results  hold  in  the  case  of  wa^es  comiiij^  in  from  +00.   We  may 


assume 


u^(i)  ^  e"^  /^  ^  H-  H^  e^  /^  ^  ,  X  —^00 


u*(x)  --  T^e"-  >^' 


.-00 


where  R^  is  the  reflection  amplitude  coefficient  and  T_^  is  the  transmission 
amplitude  coefficient  for  waves  coming  in  from  +00.  As  before,  since  (77) 
is  real,  -"i^Cx)  is  a  solution,  and  the  Wronskian  of  n,,ii^   is  constant  for  all 
values  of  x.  We  have 


u*(x)-T*  e^  /^^  ,  X  ->oo 

^L\.^ii   =-2iA  Li-i«j^: 

00 


and  therefore 


(85)  1-   |Rj^     =    |Tj^  or    |rJ2*   1t^|2  =  1. 

Another  important  relation  may  be  obtained  by  considering  the  WronskisLn 
of  u,  and  Up.   Since  u,    and  Up  are  both  so  solutions  of  (77)  we  have 


W[u^.u^  =   W[u^.u2] 


00         -00 


Then 


W[u^,u2  =  +21  /x"  T_, 
00 

Wj^u^.u^   =  f2i  ^T^ 


-00 

therefore  we  get 

(86)  '^_  =  "^^  • 
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which  states  the  reciprocity  relation  that  the  treinsmission  coefficient  for  waves 
coming  from  the  left  is  equal  to  the  transmission  coefficient  for  waves  coming 
from  the  right. 

Our  final  relation  Is  obtained  "by  considering  Wll^»^2-J*  ^®  l^a^e 


W^u^.UgJ  =  2i  y\R*  T_ 

00 


-00 

and  therefore 

(87  )  R*  T__  +  R_  T *  =  0  . 

It  is  worthwhile  noticing  that  the  functions  u.  and  u  are  unique.  Sup- 
pose there  were  two  functions  u-  ajad  wL  which  had  the  same  properties;   then 
u,-u,  would  "be  a  solution  of  (77)  which  is  outgoing  at  "both  +oo  and  -co.   For 

one-dimensional  problems  the  impossibility  of  such  non-vanishing  solutions  for 

real  positive  values  of  X  can  be  easily  established  .  Notice  also  that  u^/T.and  u^/T^ 

are  analytic  functions  of  \. 

A. The  scattering  matrix. 

We  now  consider  waves  coming  in  from  both  -oo  and  +oo.  Suppose  that 
^   (x)  contains  a  wave  a_  e  v  ^  coming  in  from  -oo  and  a  wave  a+  e~^  v^^  coming 
in  from  +oo.  Since  any  linear  combination  of  u.  and  u-  is  also  a  solution  of 
( 77  ) ,  we  have 

V(x)  =  a_U2(x)  +  a^u^(i), 
and  then  from  (81)  and  (85)  we  find  that  the  scattered  outgoing  waves  are 

^^  ^-i  A  ^  =  (^_Il_^  ^^T^)  e"^  ^  ^  ,  x->-(X, 

(88)  f- 

^_   e*^  ^^   ^  =  (a_T_+  a^R^)  e^  ^  ^  .  x^od  . 

t      * 
The  current  is,  as  stated  previously,  -r—r  W(  V  ,'/').  However  it  is  convenient 

and  physically  natural  to  think  of  the  cxirrent  as  consisting  of  incoming  and  out- 
going parts,  the  incoming  current  at  x  =+oo  for  example  being  defined  by  using  in 
the  Wronskian  not  the  complete  solution  "^  but  only  that  portion  which  is  incoming 
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at  X  =+0D  the  outgoing  CTirrent  Uq  is  proportional  to  e  '   "^  ,  the  incoming  cur- 
rent u.  is  proportional  to  e"^"/^'  ,  and  the  interlerence  terms  in  the  current 


at  X  =+00.  The  currents  may  he  interpreted  in  this  way  because  there  is  no  inter- 
ference between  the  incoming  and  outgoing  portions  of  the  current.   To  illustrate, 

ent  u 
-i/Xj 

are  proportional   to  W(u*,u^)  =  W(e^^^  ,   e^*^^^  )  and  to  W(u*,u^)  =   (e'^'^^e"^'^  ""), 
which  vanish.      Consecfuently  the  currents  at  x  =  +00  and  x  =  -00  are 

fe     W(  V  *.  ^    )       =  ^  w(a;  e^*^'  +  b*  e-^*^^  .   a^  e'^^^^  +  b     e'^^  ) 
+00  m        \  -  -  / 

+  00 

=  |-v/X      (-laj^+lbj^   ), 

^  W(  V  ^    f_)^     =  ^     W   (a:  e-^/^-    +  b:  e^^-  .   .y^-  .  b,   e-^/^0^ 
'°°     =     |A     (|a.|2-|bj2). 
It  is   seen  that  the  current  at  x  =  +00  or  x  =    -00  is  the  difference  between  the 
incoming  and  outgoing  parts,    in  the  sense  explained  above.      Moreover,   equating 
the  Wronsklans  at  +00  and  -00    , 

|/X      (|aj2+   |aj2)=|/>;   (|bj2.   |bj2    )   . 

Thus  we  have  proved  that  the  total  outgoing  current  equals  the  total  incoming 
current,  a  result  which  is  essential  for  a  consistent  interpretation  In  terms 
of  incoming  and  outgoing  currents. 

We  may  express  (88)  very  simply  in  matrix  notation.  We  have 


where  S  is  the  matrix 

(89)  /t^  R_ 

R^.  T_ 
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S  is  called  the  collision  or  scattering  matrix  and  is  defined  so  as  to  equal  the 
unit  matrix  when  q  is  zero  in  (77  )•   From  (83),  (85),  ancl  (87)  we  obtain  the 

relation  . 

S  S  =  I, 

where  S     is   the  complex  conjugate  of  the  transpose  of  S,    and  I   is  the  unit  matrix. 
This  relation  shows   that  S  is  unitary.     This  property  of  S  coxild  also  have  teen 
deduced  from  the  fact  that 

1^1^-  \K\^=  i^j'*  \^y ' 

Note  that  from  (86)  and  (87)  it  follows  that 

lTj2=  iTj^a^d  lRj2=  lHj2  , 
which  shows  that  the  reflection  and  transmission  coefficients  do  not  depend  on 
whether  the  wave  comes  from  the  left  or  right. 

Up  to  this  point  we  have  assumed  that  X  is  real  in  equation  (77)  and 
the  subsequent  ansdysis.  If  X  is  complex,  u^  eind  Up  are  no  longer  solutions  of 
the  equation  (77)  ^or  ^^®  same  eigenvalue  X  as  are  tL.  and  u_;   consequently  the 
relations  (83), (85),  and  (87)  are  no  longer  valid.  Only  the  relation  (86)  remains. 

We  may,  however,  continue  to  define  the  matrix  S  of  (89)  formally  as 
the  analytic  continuation  into  the  complex  X-plane  of  the  scattering  matrix  ob- 
tained for  real  values  of  X.  Because  the  relations  (83),  (85),  and  (87)  are  not 
valid,  S  will  no  longer  be  unitary  for  complex  values  of  X.   It  can  easily  be 
shown  that  T  ,  T  ,  R  ,  and  R  are  ginalytic  functions  of  X  in  the  upper  and  lower 
half-planes.   Thus  the  only  possible  singularities  of  the  scattering  matrix  sire 
on  the  real  axis. 

We  shall  show  that  the  singularities  of  the  scattering  matrix  contain 
the  singularities  of  the  Green's  function  G(x,x',X)  for  the  operator  L  and 
therefore  also  for  the  eigenvalues  ol  L.  The  solution  of  (77)  ^ich  for  com- 
plex X  is  quadratically  integrable  at  x  =  +00  is  proportional  to  Up(x,X),  i.e. 
to  the  solution  which  in  this  section  has  been  called  outgoing  at  +00,   In  fact 
the  function  Up/T  behaves  like  e  ^   at  +00.   Similarly  the  function  u.^/T^  behaves 


-i/Xx 


like  e  "'  at   -00,     Thus^  since 


W[u^/T^.  U2/T_]   =  ^     ^\3^.^2^        =  ^ 


00 
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we  have  as   in  section  VII 

-i/Xx       i/Xi'                    l/X  (x'-x) 
G(x,x'.\)  =  -     -±- 5 =  -     -? T  U)  , 

WL^/T+.V^-^  21  /\ 

for  X   <  X*,   X  approaching  -oo  emd  x'   approaching  +oo   . 

We  shall   now  derive  integral  expressions   for  the  elements  of  the 
scattering  matrix.      Consider  the  function  Up.      We  know  that 

where  0p(x)  is   outgoing  at  both  +oo  smd  -oo   .      This  expression  is   substituted  for 
Uo  in  (77)  and  we  get 

-     p—     -     X^p     =  -qu     . 

Now  we  use  the  outgoing  wave  Green's  function  Eq.    (57)  for  the  left  side  and 
obtain 

2         21/5:      7  ^ 

or 

00 

-  iA^ -.     1     A      i/>r|x-x'|    ,  ,.   ,  ,v^  , 

Uj  =  e  *         +     /  e   »      '  '    q(x«)u  (x')dx»    . 

'^  2i/\    »/-oo  '^ 

Now  let  X  approach  first  -oo  and  then  +oo.  We  find  that 

i/X  X   .   „       -iA^-     i/X3c.      g-iAi    Z'        i/^^, 
e"         +Re'         =e''  +     /       e   '         ou„dx' 

2iyX     ^-00  ^ 

and 


iy/X 


00 


^    i/xx_  lyxx  .  e^y^  r    -i/xi'      ^  , 

T     e   "         =  e  *  +  I        Q     ^  qu     dx'; 

2i  y\      ^-00  2 

consequently 


T     =  1  +  /        e      *  qup   dx 


2i/X     ^''-00 


1  /•        i,/Xx  . 
/       e   *         qu-  dx     . 

2iyX     ^oo  '^ 
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Similarly  we  obtain 

00      rr 

T  =  1  +  -^^—  /    e^^^  ^  qu-  dx 

*  2iyx  •'-oo  ^ 

(91) 


1 r'^    -lyxx 

E^  =  /    e  *     qiL,  dx 


2i  /X   ^-00 


Reference  13  gives  a  much  expanded  treatment  of  the  scattering  prohlem. 

In  this  work,  it  is  shovra  that  it  is  also  possible  to  derive  results  such  as  (90) 

13  Ik 
and  (91)  hy  means  of  formal  symbolic  operator  techniques      ,  We  have  pre- 
ferred, however,  to  demonstrate  (90 )  ajid  (9I)  directly  from  the  differential 

equation  in  order  to  avoid  any  difficulties  about  the  meajiing  of  the  operators 

13  lU 
that  must  otherwise  be  symbolically  defined  ■^'   .  Ve  intend,  in  our  next  re- 
port, to  consider  the  scattering  matrix  for  a  three-dimensional  particle,  in 
which  case  the  corresponding  one-dimensional  (radial)  equation  covers  the  semi- 
infinite  interval  0  £  r  <oo  instead  of  the  fully  infinite  interval  which  we 


have  been  considering.   There  exist  numerous  papers  on  the  subject  of  the 
scattering  matrix  in  the  three-dimensional  case;   thos 
contain  results  and  references  which  are  of  interest. 
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scattering  matrix  in  the  three-dimensional  case;   those  by  Ma   and  van  Kampen 
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